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Abstract. Resonant X-ray emission spectroscopy (XES) measurements at Fe L2.3- 
edges and electronic structure calculations of LiFeAs and NaFeAs are presented. 
Experiment and theory show that in the vicinity of the Fermi energy, the density 
of states is dominated by contributions from Fe 5(i-states. The comparison of Fe ^^,3 
XES with spectra from related FeAs-compounds reveal similar trends in energy and 
the ratio of the intensity of the L2 and Lg peaks {1{L2) /l{Lg) ratio). The 1{L2) /liLg) 
ratio of all FeAs-based superconductors is found to be closer to that of metallic Fe than 
that of the strongly-correlated FeO. We conclude that iron-based superconductors are 
weakly- or at most moderately-correlated systems. 
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1. Introduction 

Recently, a new class of iron-based superconductors was discovered. LiFeAs and NaFeAs 
were found to be superconducting with Tc = 18K and 9K, respectively [1-4]. These 
compounds have an unexpected lack in magnetic order at all temperatures whereas 
in i?EOFeAs (where RE is a rare earth element) and >l£'Fe2As2 (where AE is an 
alkaline earth element) the Fe magnetic moments adopt a coUinear antiferromagnetic 
(c-AFM) order at low temperatures [2]. In contrast with other iron-based compounds, 
no further doping is necessary to induce superconductivity and the spin-density wave 
(SDW) state appears to be notably absent from these new systems [5]. This suggests 
LiFeAs and NaFeAs can offer important insights to understanding the mechanism of 
superconductivity in iron-based superconductors. The absence of SDW transitions 
and the relatively low Tc in comparison with i?£'OFeAs and 74EFe2As2 iron-based 
superconductors make these two compounds possible candidates for conventional BCS 
superconductors [6] . However the estimation of the electron-phonon coupling parameter 
A from band structure calculations [6] gives low values of A = 0.29 and 0.27 for 
LiFeAs and NaFeAs respectively, which are too weak to account for conventional BCS 
superconductivity in this class of superconductors. Experimental studies of electronic 
structure of LiFeAs and NaFeAs superconductors and their comparison with that of 
i^E'OFcAs and y4i?Fc2As2 iron-based superconductors are important in resolving this 
puzzle. In this manuscript we present and analyze X-ray emission spectroscopy (XES) 
measurements at the Fe L^^^-edge of LiFeAs and NaFeAs. These measurements probe 
the occupied Fe 3d density of states (DOS). We compare the measured spectra with our 
electronic structure calculations of the valence structure. We contrast our findings with 
our previous studies of i2£'0FeAs {RE = La, Sm) [7] and CaFe2As2 [8]. 

2. Experimental and Calculation Details 

LiFeAs and NaFeAs were synthesized by the reaction of stoichiometric quantities of 
elemental Li (or Na), Fe and As. Fe and As powders were ground together and added to 
pieces of Li (or Na) in a tantalum tube, which was then sealed by welding under 1 atm 
argon gas; the mixture was heated to 800 °C for 2 days. Prehminary characterization 
of the resulting product by powder X-ray diffraction revealed pure LiFeAs and NaFeAs 
phases. For details of sample preparation see [1] and [4]. 

The X-ray emission measurements of LiFeAs and NaFeAs were performed at the 
soft X-ray fluorescence endstation at Beamline 8.0.1 of the Advanced Light Source 
at Lawrence Berkeley National Laboratory [9]. The endstation uses a Rowland circle 
geometry X-ray spectrometer with spherical gratings and an area sensitive multichannel 
detector. We have measured resonant and non- resonant Fe L2^3 {3d4s — > 2p transition) 
X-ray emission spectra (XES). The instrumental resolving power (E/AE) for Fe L2,3 
spectra was 10^. All spectra were normalized to the incident photon current using a 
highly transparent gold mesh in front of the sample to correct for intensity fluctuations 
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Compound Fe-Fe NN(Fe) Fe-R NN(R) Pt. Grp. 



LaOFeAs [12] 


2.85 


4 


2.41 


4 


-42m 


CaFeaAsa [13] 


2.74 


4 


2.32 


4 


-42m 


NaFeAs 


2.79 


4 


2.44 


4 


-42m 


LiFeAs 


2.67 


4 


2.41 


4 


-42m 


Fe (bcc) [14] 


2.54 


8 






m3m 


FeO [15] 


3.06 


12 


2.16 


6 


m3m 



Table 1. Structural parameters for various Fe compounds. Here "Fe-Fe" refers to the 
distance between iron atoms, "Fe-R" refers to the distance between the iron atom and 
the anion, both in A. NN() refers to the number of nearest neighbours of the specified 
type. The final column, "Pt. Grp." refers to the point group around the Fc atom. 

in the photon beam. The excitations for the XES measurements were determined from 
Fe 2p X-ray absorption spectroscopy (XAS) measurements (using total electron yield 
mode); the chosen energies corresponded to the location of the L2 and L3 thresholds, 
an energy between them, and an energy well above the threshold. 

All density of states calculations were performed within the full-potential linear 
augmented plane- wave method as implemented in WIEN2k code [10]. For the exchange- 
correlation potential we used the generalized gradient approximation in the Perdew- 
Burke-Ernzerhof variant [11]. The Brillouin zone integrations were performed with a 12 
X 12 X 7 special A;-point grid for LiFeAs and a 12 x 12 x 6 special /^-point grid for 
NaFeAs. and RJ^^^Kmax = 7 (the product of the smallest of the atomic sphere radii Rmt 
and the plane wave cutoff parameter Kmax) was used for the expansion of the basis set. 
The experimental values of the high-temperature lattice constants and atomic positions 
are used [1,4]. Both compounds have a tetragonal crystal lattice in the P4/nmm space 
group. In NaFeAs, the lattice constants are a = 3.9494 A and c = 7.0396 A, whereas 
in LiFeAs, a = 3.7754 A and c = 6.3534 A (for a summary of structural parameters see 
Table [H). The atomic sphere radii were chosen as Rnq = 2.5, Rpe = 2.3, Ras =2 .04 

a. u. and Rn = 2.36, Rpe = 2.28, Ras = 2.02 a.u. for NaFeAs and LiFeAs respectively. 
The sphere radii were chosen in such a way that the spheres were nearly touching. 

3. Results and Discussion 

The measured XES and XAS spectra are shown in Figure [H The LiFeAs and NaFeAs Fe 
Lg 5 XES (Fig. [1], bottom panels) indicate two main bands located around 705 and 718 
eV, these correspond to the Fe Lg {SdAs 2^3/2 transitions) and Fe Lg {SdAs — > 2pi/2 
transitions) normal emission lines separated by the spin-orbital splitting of Fe 2p-sta.tes. 
The non-resonant Fe Lg^^ XES (XES curve a) for both materials lacks the low-energy 
satellite structure typical for correlated systems (for instance for FeO [16]) and the main 
peak is sharp and similar to metallic iron [17]. The resonant XES spectra (XES curves 

b, c, d in the bottom panels of Figure [1]) show no energy-loss features; this indicates 
that even resonant Fe Lg^^ XES probes mainly the partial DOS in these materials. The 
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Figure 1. Summary of spectra for LiFeAs (left side) and NaFeAs (right side). The 
excitation energies for resonant Fe ^2,3 XES are indicated by arrows in the XAS 
spectra in the left panels. 



XAS spectra (top panels of Figure [T]) indicate the absorption thresholds for Fe 2pi/2 
and 2p3/2 electrons, and are used to determine the appropriate excitation energies for 
resonant XES. The detailed features in XAS spectra reveal information about the Fe 3d 
unoccupied states (the conduction band), however due to the greater effective atomic 
potential when a 2p core electron is removed the states in a XAS spectra are greatly 
distorted from those in the unperturbed crystal; this is referred to as the "core-hole 
effect". 

The integral of the L2 and Lg peaks in a non-resonant XES measurement are related 
to the population of the 2pi/2 and 2p3/2 states respectively. For free atoms the ratio 
of the integrals of the and L3 XES peaks (the I(i/a)/I(L5) ratio) should be equal 
to |. In metals the radiationless L2L3M45 Coster-Kronig (C-K) transitions strongly 
reduce the l{L2)/l{L3) ratio [18]. Figure [2] illustrates that the I(L2)/I(Z5) ratio is 
almost identical for all previously mentioned FeAs-compounds. This ratio is closer to 
that of metallic Fe than to that of FeO, and indicates that the Fe 3d states in FeAs- 
compounds are much less correlated then those in FeO. Further, the low energy edge 
of the L3 peaks in the FeAs-compounds and metallic Fe lack the prominent satellite 
that is present in FeO. As previously mentioned this also supports the conclusion that 
the Fe 3d electrons in FeAs-compounds are largely itinerant. For metallic compounds, 
the Fe 3d bandwidth should decrease with increasing Fe — Fe distance and decreasing 
number of Fe — Fe neighbours (see Table [1]), and this is somewhat demonstrated in 
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Figure 2. Comparison of tiie L2.3 XES spectra for metallic Fe, LiFeAs, LaOFeAs [7], 
CaFe2As2 [8], NaFeAs, and FeO. Note the low energy shoulder s appears in FeO but 
not in any of the other spectra. The left inset shows the full-width at half maximum 
(FWHM) of the Fe L3 peak. The right inset shows the I(Lg)/I(L3) ratios for metaUic 
Fe, LaOFeAs, CaFe2As2, and correlated FeO. The 1{L2)/I{L3) ratio was calculated 
from the ratios of the integral under the L2 and L3 peaks, respectively. 



the full-width at half maximum (FWHM) of the Fe Ls peak (Figure [2], top left inset), 
which reduces from ~3.3 eV for Fe metal to ~3.0 eV for NaFeAs. Except for LaOFeAs, 
which we would expect to have a smaller FWHM than NaFeAs, the trend in decreasing 
FWHM qualitatively matches the trend in increasing Fe — Fe distance in Table [1] for 
the metallic compounds. It should be noted that these changes in FWHM are less then 
the instrumental resolution by roughly a factor of two and that core-level spectroscopy 
is not an appropriate technique to probe this effect - so we do not attempt to draw any 
conclusions from this trend. There is, however, a clear difference in the FWHM between 
the FeAs-compounds and FeO. 

The calculated electronic structure of LiFeAs and NaFeAs are presented in Figure 
[31 Our calculations match those performed in Reference [6]. The calculated partial 
density of states (DOS) distribution is quite similar for both compounds. In particular, 
the density of states in the vicinity of the Fermi level (0 to -2 eV) is dominated by 
contributions from Fe Sd-states (region a) and As ^j)-states mixed with Fe 3d,4s-states 
(region b) are located at the bottom of the valence band (-2 to -5.5 eV). The alkali earth 
2s, ^s-states (for Li and Na, respectively) provide a minimal but consistent contribution 
from to -5.5 eV, these states hybridize with the Fe 3d-sta.tes. 
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Figure 3. Calculated DOS for LiFeAs and NaFeAs. The dotted lines in the Fe 3s, 3p 
DOS plot refer to the 3s states increased by a factor of 10. The As 4^ states are 
separated from the 4P states so they are plotted with the same line style. The y-axis 
in the total DOS plot is in units of [states/eV/unit cell]. Two regions, a and b are 
identified in the total DOS plot. Fe 5d-states dominate in region a. Roughly even 
contribution from Fe 3d-states and As ^p-states occurs in region b. 



To compare the XES spectra with the calculated DOS in a meaningful way, we 
separated the L3 band in the XES spectra into pseudo-Voigt components (of the 
form given in [1]). As discussed, XES probes the partial occupied density of states; 
however the DOS features are broadened by both the instrumental resolution (Gaussian 
in nature) and the core- hole lifetime [19] (Lorentzian in nature). We only fit the 
L3 band since the L2 band is basically the same partial occupied density of states 
with poorer statistics. The position fii, amplitude Ai, and Lorentz broadening Fj for 
each component are determined by least-squares fitting. Here we set the Gaussian 
a parameter to the position /ij divided by the instrumental resolving power (E/AE), 
and kept a consistent mixing factor 77 for all component peaks. We calculated best- 
fit curves with one to six pseudo-Voigt components and the quality of fit parameter 
F" = {fdata{x) — f fit{x))^ was examined to determine the "simplest best fit". 

Since least-squares fitting requires initial estimates for the fitted parameters, fitting 
was conducted several times for a given number of pseudo-Voigt components with range 
of different initial positions /Xj; this quantity is the hardest to fit accurately and is also 
the one we are most interested in. The initial estimates for A,, F,, and r] were set to 
1.0, 1.0, and 0.5 respectively in all cases, these quantities are fitted quite well regardless 
of the initial estimate. We obtained consistent fits for several different initial estimates 
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(1) 



For LiFeAs four or more pseudo-Voigt components produced fits with F" values 
within 10% of each other while the F" for the best fit with 3 components was ~25% 
greater than the 4 component fit (see Figure HI top panel, inset). Likewise, 3 or more 
pseudo-Voigt components produced a consistently good fit for NaFeAs, whereas the 
F" for the best 2 component fit was ~50% greater than that of the 3 component fit. 
The fit for NaFeAs was not as good as for LiFeAs because the NaFeAs data had more 
noise. The fitted curve for LiFeAs matches the measured Fe spectrum, and the 
components are located around the main features in the Fe 3d DOS calculation (Figure 
m top panel). Only the main Fe 3d DOS feature at 705.5 eV is sharp and isolated 
enough to potentially provide a pure pseudo-Voigt spectrum, so we expect the fitted 
components at lower energies to have higher amplitudes relative to the pseudo-Voigt at 
705.5 eV than a comparison of the representative heights of the Fe 3d DOS features at 
the Fermi energy and elsewhere in the valence band would indicate. For example, the 
ratio between the heights of the Fe 3d DOS feature plotted at 705.5 eV and the feature 
plotted at 702.8 eV in the top panel of figure H] is much greater then the ratio between 
heights of the pseudo-Voigt component at 705.5 eV and the pseudo-Voigt component at 
703.1 eV (see Table [2]), which is not unexpected since there are many other Fe 3d DOS 
features at ~703 eV within the instrumental resolution of each other, and all of them 
would contribute to the spectrum, but would not be independently resolvable. 

The fitted curve for NaFeAs matches the measured Fe L3 spectrum, but only one 
of the pseudo-Voigt components is in good agreement with a Fe 3d DOS feature (see 
Figure HI centre panel) . It is likely that the poorer quality of the NaFeAs XES spectrum 
compared to the LiFeAs XES spectrum is to blame here; it is possible to obtain a good 
fit using 4 pseudo-Voigt components fixed at 705.5 eV, 704.8 eV, 703.0 eV and 701.2 eV 
as the features in the Fe 3d DOS suggest, however our aim was to identify features in 
the DOS from the XES spectra without requiring input from calculations. 

Resonant XES spectra usually have better signal-to-noise ratio, and since these 
materials do not show significant energy-loss features in the resonantly excited spectra, 
we can attempt to fit the resonantly excited Fe L3 spectrum by the method outlined 
above. In this case, 4 pseudo-Voigt components provide the optimal fit (Figure HI inset 
in bottom panel) and these components give a much better agreement to the Fe 3d DOS 
than the 3-component non- resonant XES fitting did (see Figure HI bottom panel). The 
fitted curve matches the measured spectrum except at the edge above 708 eV. This is 
not unexpected, since 708 eV is the excitation energy for this spectrum (see Figure [H 
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Figure 4. Comparison of calculated DOS and measured spectra. The top panel shows 
the L3 peak of the Fe L2.3 spectrum, the Fe 3d DOS, best-fit curve, and the 4 pseudo- 
Voigt components contributing to the best-fit curve for LiFeAs. The centre panel shows 
the Lg peak of the Fe L2.3 spectrum, the Fe 3d DOS, best-fit curve, and the 3 pseudo- 
Voigt components contributing to the best-fit curve for NaFeAs. The bottom panel 
shows the resonantly excited Fe L3 spectrum, the Fe 3d DOS, best-fit curve, and the 
3 pseudo-Voigt components contributing to the best-fit curve for NaFeAs. The insets 
show the normalized fit parameter F" for fits with different numbers of pseudo-Voigt 
components for each material. Note that the maximum F" has been scaled to 1.0 in 
each case - for example in LiFeAs the worst fit (with 1 pseudo-Voigt function) falls far 
outside the range of the plot. The estimated Fermi level is indicated in each plot. 



top right panel) and the elastically scattered X-ray peak will distort the spectrum near 
708 eV - in particular by increasing the amplitude of the feature with less broadening 
(the elastic scatter has a purely Gaussian profile) then a pseudo-Voigt with parameters 
consistent with the remainder of the spectrum will be able to satisfy. The Fe L3 XES 
band from the other measured spectra (curves b and c in Figure [1] for both LiFeAs 
and NaFeAs) produced results consistent with those discussed above. In general, non- 
resonant XES should provide the best results if the signal-to-noise ratio is high enough, 
XES resonant with the L3 feature may be used if there are no energy-loss features or 
prominent scattering features. 

The Fermi levels in FigurelHwere estimated by aligning the calculated DOS with the 
fitted pseudo-Voigt components, so the agreement between DOS and the first pseudo- 
Voigt component is manufactured. One method of estimating the Fermi level from 
emission spectra is to use the peak of the second derivative [20]. We have estimated 



Identifying valence structure in LiFeAs and NaFeAs with core-level spectroscopy 9 





Vi 


Ai 


fli [eVJ 


ij [eVJ 


Lii^eAsi 


(J.47 


1 AO 

1.48 


705.5 


1 OA 

1.34 


T *T7 A 

LireAs2 




1 A A 

1.44 


704.4 


9.1/ 


T * 1~i A 

Lii'eAss 




0.96 


703.1 


5.24 


LiFeAs4 




0.29 


701.3 


4.15 


AT TTi A TV 

NaFeAsf 


0.37 


2.07 


705.5 


2.27 


IN Ac ei\k52 




1.410 




^ 77 


NaFeAsf 




0.53 


702.5 


4.19 


NaFeAsf 


0.39 


1.28 


705.8 


1.25 


NaFeAsf 




1.05 


704.9 


1.61 


NaFeAsf 




1.12 


703.5 


3.31 


NaFeAsf 




0.43 


701.7 


4.62 



Table 2. Fit results for the 4 pseudo-Voigt peaks for fitting the Fe XES of LiFeAs 
and the 3 and 4 pseudo-Voigt peaks for fitting the nonresonant and resonant Fe L3 XES 
of NaFeAs (NaFeAsf^ and NaFeAsf respectively). Note there is only one mixing- factor 
?7 for all components in the same spectrum. 

Fermi levels of 705.9 eV, 706.0 eV, and 706.3 eV from aligning DOS and fit components 
and 706.5 eV, 706.8 eV, and 706.7 eV from the peak in the second derivative of the 
appropriate Fe L3 XES for for LiFeNa, non-resonant NaFeAs, and resonant NaFeAs, 
respectively. Since the main peak at 705.5 eV is so close to the Fermi level, the Fe 
3d occupied states suffer an abrupt cut-off rather than a gradual decline at the Fermi 
level, and the Ls XES portions at higher energies than that are due only to spectral 
broadening. It is therefore not unexpected that the XES second-derivative estimates of 
the Fermi level are greater than the curve-fit alignment Fermi level estimate by a shift 
of roughly the instrumental energy resolution. 

To summarize, we have studied the electronic structure of LiFeAs and NaFeAs via 
resonant and non-resonant XES spectra and DFT calculations using the generalized 
gradient approximation. We demonstrate that pseudo-Voigt curve fitting of XES 
measurements without knowledge of the electronic structure can give good agreement 
with prominent features in the calculated valence band. The results from DFT 
calculations and the comparison of Fe Ls FWHM and \{L2) /\{Ls) ratio with standard 
reference compounds suggests the Fe 3d states are mostly itinerant in nature. The 
comparison of our results with previous studies on LaOFeAs and CaFe2As2 show that 
the main features of electronic structure found here are general for all studied FeAs- 
systems: Fe 3d-states dominate near the Fermi level and As 4p-states are concentrated 
at the bottom of the valence band. We conclude that all FeAs-systems studied herein 
are weakly or moderately correlated systems. 
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